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We present the first measurement of multiplicity and pseudorapidity distribution of pliotons in 
tlie pseudorapidity region 2.3 < rj < 3.7 for different centralities in Au + An collisions at ^/sNN 
= 62.4 GeV. We find that the photon yield in this pseudorapidity range scales with the number 
of participating nucleons at all collision centralities studied. The pseudorapidity distribution of 
photons, dominated by neutral pion decays, has been compared to those of identified charged pious, 
photons, and inclusive charged particles from heavy ion and nucleon-nucleon collisions at various 
energies. The photon production in the measured pseudorapidity region has been shown to be 
consistent with the energy and centrality independent limiting fragmentation scenario. 

PACS numbers: 25.75.Dw 



One of the primary goals of the heavy-ion pro- 
gram at the Relativistic Heavy-Ion Collider (RHIC) 
at Brookhaven National Laboratory is to search for 
the possible formation of Quark-Gluon Plasma 0. 
Important information about the dynamics of par- 
ticle production and the evolution of the system 
formed in the collision can be obtained from various 
global observables, such as the multiplicity of pho- 
tons and charged particles. At RHIC energies, the 
particle production mechanisms could be different 
in different regions of pseudorapidity (77) |^ . At 
midrapidity a significant increase in charged particle 
production normalized to the number of participat- 
ing nucleons {Npart) has been observed for central 
Au+Au collisions compared to peripheral Au+Au 
and p+p collisions 4]. This has been attributed to 
the onset of hard scattering processes, which scale 
with the number of binary collisions. Alternatively, 
in the Color Glass Condensate Q picture of par- 
ticle production at midrapidity, the centrality de- 
pendence could reflect increasing gluon density due 
to the decrease in the effective strong coupling con- 
stant. However, the total charged particle multiplic- 
ity per participant pair, integrated over the whole 
pseudorapidity range, is found to be independent of 
centrality in Au-t-Au collisions 0- 

It is also observed that the number of charged par- 
ticles produced per participant pair as a function of 
V - Ybeam, whcrc Ybeam IS the beam rapidity, is in- 
dependent of beam energy j2| • This phenomenon is 
known as limiting fragmentation (LF) There 
have been contradictory results reported from in- 
clusive charged particle measurements regarding the 
centrality dependence of the LF behavior, results 
from PHOBOS show a centrality dependence 0, 
while those from BRAHMS show a centrality inde- 
pendent behaviour j2i]. The centrality dependence 
at forward rapidities has been attributed to nuclear 
remnants and baryon stopping. The role of a new 
mechanism of baryon production 7] also needs to be 
understood. Further insight into this question can 
be obtained by studying the centrality, beam energy 
and system size dependence of LF phenomena with 
identified particles. Beam energy independence of 



LF for identified pions has been found in e+e^ col- 
lisions 3 . 

Photons are produced in all stages of the system 
created in heavy ion collisions. They do not inter- 
act strongly with the medium and carry information 
about the history of the collision. Since inclusive 
photon production is dominated by photons from 
the decay of tt'^'s, measurement of the multiplicity of 
photons is complementary to the charged pion mea- 
surements. The forward rapidity region in heavy 
ion collisions, where the present measurements have 
been carried out, constitutes an environment that 
precludes the use of a calorimeter due to the high 
level of overlap of fully developed showers. The only 
measurements of photon multiplicity distribution in 
the forward rapidity region reported to date are from 
a preshower detector 9] at the SPS, resulting in the 
study of various aspects of the reaction mechanism 
in heavy ion collisions 0, ^| . 

In this Letter we present the first measurement of 
photon production at the forward rapidities (2.3 < 
77 < 3.7), carried out by the STAR experiment using 
a highly granular preshower photon multiplicity de- 
tector (PMD) 12] in Au -I- Au collisions at y^J/vw = 
62.4 GeV. The STAR experiment consists of several 
detectors to measure hadronic and electromagnetic 
observables The minimum bias trigger is ob- 

tained using the charged particle hits from an array 
of scintillator slats arranged in a barrel called the 
Central Trigger Barrel surrounding the Time Projec- 
tion Chamber (TPC) and two zero degree hadronic 
calorimeters at ± 18 m from the detector center 
A total of 334000 minimum bias events, correspond- 
ing to to 80% of the Au+Au hadronic interaction 
cross section, have been selected with a collision ver- 
tex position of less than 30 cm from the center of the 
TPC along the beam axis. The centrality determina- 
tion in this analysis uses the multiplicity of charged 
particles in the pseudorapidity region | 77 | < 0.5, as 
measured by the TPC H^. 

The PMD is located 5.4 meters away from the cen- 
ter of the TPC (the nominal collision point) along 
the beam axis. It consists of two planes (charged 
particle veto and preshower) of an array of cellular 
gas proportional counters [13I ■ A lead plate of 3 ra- 
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diation length thickness was placed between the two 
planes and was used as a photon converter. The sen- 
sitive medium is a gas mixture of Ar and CO2 in the 
ratio of 70%:30% by weight. There are 41472 cells 
in each plane, placed inside 12 high voltage insu- 
lated and gas-tight chambers called super modules 
(SM). A photon traversing the converter produces 
an electromagnetic shower in the preshower plane, 
leading to a larger signal spread over several cells as 
compared to a charged particle which is essentially 
confined to one cell [IjI- In the present analysis, 
only the data from the preshower plane have been 
used. 

The cell-wise response is obtained by using the 
ADC distributions of isolated cells. The ADC dis- 
tribution of an isolated cell may be treated as the 
response of the cell to charged particles For 
most of the cells this response followed a Landau 
distribution. We used the mean of the ADC distri- 
bution of isolated cells to estimate and correct the 
relative gains of all cells within each SM. The cell- 
to-cell gain variation within a SM varied between 10 
- 25% for different SMs. 

The extraction of photon multiplicity proceeds in 
two steps involving clustering of hits and photon- 
hadron discrimination. Hit clusters consist of con- 
tiguous cells. Photons are separated from charged 
particles using the following conditions : (a) The 
number of cells in a cluster is > 1 and (b) the clus- 
ter signal is larger than 3 times the average response 
of all isolated cells in a SM. The choice of the above 
condition is based on a detailed study of the detec- 
tor response using Monte Carlo simulations. The 
number of selected clusters, called 7 — like clusters 
(Nj-iike), in different SMs for the same rj coverage is 
used to evaluate the effect of possible non-uniformity 
in the response of the detector. 

To estimate the number of photons (N-y) from the 
detected N-y-uke clusters we evaluate the photon re- 
construction efficiency (e-y) and purity (/p) of the 
7 — like sample defined [lO| as = Njd^' /N^ and 
fp = N2{1^ /N^-iike respectively. iVj,'J^ is the num- 
ber of photon clusters after applying the photon- 
hadron discrimination conditions. Both and fp 
are obtained from a detailed Monte Carlo simu- 
lation using the HIJING event generator (version 
1.382) [3 

with default parameter settings and the 
detector simulation package GEANT fl7l| , which in- 
corporates the full STAR detector framework for the 
period this data was taken. The lower limit of pho- 
ton pt acceptance in the PMD is estimated to be 20 
MeV/c. Both e-y and fp vary with pseudorapidity 
and centrality. This is due to variations in parti- 
cle density, upstream conversions and detector re- 
lated effects. The highest occupancy is about 12% 
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FIG. 1: Minimum bias A^^ distribution. Comparison 
with HIJING and AMPT models are shown. Horizontal 
bars indicate the errors. The photon multiplicity distri- 
bution for top 5% central events is shown in open circles. 
The solid curve is the fit by a Gaussian. 



and maximum percentage of split cluster is 9%. The 
photon reconstruction efficiency is determined from 
simulations to increase from 42% to 56% in central 
collisions and from 42% to 70% in peripheral col- 
lisions as rj increases from 2.3 to 3.7. The purity 
of the photon sample ranges from 55% to 62%, and 
from 63% to 70% for central and peripheral colli- 
sions respectively as we increase rj within the above 
range. 

The systematic errors on the photon multiplicity 
(iV^) are due to (a) uncertainty in estimates of 
and fp values, arising from splitting of clusters and 
the choice of photon-hadron discrimination condi- 
tions and (b) uncertainty in iV^ arising from the 
non-uniformity of the detector primarily due to cell- 
to-cell gain variation. The error in N.^ due to (a) is 
estimated from Monte Carlo simulations to be 9.8% 
and 7.7 % in central and peripheral collisions respec- 
tively. The error in due to (b) is estimated using 
average gains for normalization and by studying the 
azimuthal dependence of photon density of the de- 
tector in a 77 window. This is found to be 13.5% for 
central and 15% for peripheral collisions. The to- 
tal systematic error in is ~ 17% for both central 
and peripheral collisions. The errors are obtained by 
adding systematic and statistical errors in quadra- 
ture and are shown in all the figures. The statistical 
errors are small and within the symbol sizes. 

Fig. n shows the minimum bias distribution of 
N-^ along with results from HIJING events passed 
through detector response (henceforth referred to as 
HIJING) and AMPT ^ models. The sharp drop 
in HIJING results at higher is due to lack of 
statistics. The HIJING model is based on perturba- 
tive QCD processes which lead to multiple jet pro- 
duction and jet interactions in matter. The AMPT 
model is a multi-phase transport model which in- 
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FIG. 3: Variation of A^.^ per participant pair in PMD 
coverage (2.3 < r) < 3.7) as a function of Npart- The 
lower band reflects uncertainties in Npart calculations. 



eludes both initial partonic and final hadronic in- 
teraetions. We observe that HIJING underpredicts 
the measured photon multiplicity whereas AMPT 
slightly overpredicts the total measured photon mul- 
tiplicity for central collisions. Within the errors, the 
two models are in agreement with the measurement. 
The top 5% central photon multiplicity distribution 
(open circles) is fitted by a Gaussian with a mean of 
252. 

Fig. 13 shows the pseudorapidity distribution of 
photons for various event centrality classes. The 
results from HIJING are systematically lower com- 
pared to data for mid-central and peripheral events. 
The results from AMPT compare well with the data. 

Fig-Elshows the variation of total number of pho- 
tons per participant pair in the PMD coverage as a 
function of the number of participants. Npart is ob- 
tained from Glauber calculations [l^. Higher val- 
ues of Npart corresponds to central collisions. We 
observe that the total number of photons per par- 
ticipant pair is approximately constant with central- 
ity. The values from HIJING are lower compared to 
the data. The values from AMPT agree fairly well 
with those obtained from the data. Approximate lin- 
ear scaling of N^ with Npart in the rj range studied 
indicates that photon production is consistent with 
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FIG. 4: (Color Online) (a) Variation of 
to Npart with rj 



dN-, 
drf 



normalized 



ybeam for different collision energy and 
centrality. Also shown for pp collisions, (b) same as 
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nucleus-nucleus collisions being a superposition of 
nucleon-nucleon collisions. 

In Fig. 0] we present the energy and centrality de- 
pendence of LF for inclusive photons and charged 
particles. Fig. EJa) compares the photon pseudora- 
pidity distributions for central (0-5%) and peripheral 
(40-50%) Au + Au collisions at .JsWn = 62.4 GeV, 
with the top SPS ener gy c entral (0-5%) photon data 
for Pb -I- Pb colhsions [13 as a function of ?7 - ybeam- 
Also shown are the photon data from pp collisions 
at ./FWn = 540 GeV 0. In Fig. E^b) we show 
the charged particle pseudorapidity distributions for 
central (0-6%), peripheral (35-40%) Au + Au colli- 
sions at ^Jsnn = 200 GeV and central data at ^Jsjsin 
= 130 GeV from PHOBOS and BRAHMS H as 
a function of ry-ybeom- Also shown are the charged 
particle data from pp and pp collisions at ^/snn = 53 
and 200 GeV [13 • We observe in Fig.g^a) that pho- 
ton results from the SPS and RHIC are consistent 
with each other, suggesting that photon production 
follows an energy independent LF behavior. Similar 
energy independent LF behavior had been observed 
for charged particles 0, Q . This can again be seen 
in Fig.^b) from the comparison of charged particle 
-q distributions from PHOBOS for y^sJm = 130 and 
200 GeV and BRAHMS at = 130 GeV . 

In Fig. El^a) we also observe that within the mea- 
sured 7y range the photon distribution as a function 
of V - Ybeam is independent of centrality. However, 
in Fig. 01b) it is observed that the charged particle 
distribution as a function of rj - ybeam is dependent 
on centrality 0] . This centrality dependent behavior 
of LF observed by PHOBOS is most prominent at 
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FIG. 5: (Color Online) Estimated from nor- 

malized to Npart, as compared to '^"'^^^ normalized to 
Npart, as a function of y - jheam for central collisions at 
various collision energies. 



the lower energy of = 19.6 GeV Q. The cen- 

trality dependence has been speculated to be due to 
nuclear remnants and baryon stopping 0, 01 ■ The 
dependence of LF on the collision system is most 
clearly seen in the comparison between results from 
heavy ion collisions with those from pp and pp col- 
lisions. We observe in Fig. EJa) that the photon 
results in the forward rapidity region from pp col- 
lisions at ^/SNN — 540 GeV are in close agreement 
with the measured photon yield in Au-|-Au collisions 
at y^STvw = 62.4 GeV. However the pp and pp inclu- 
sive charged particle results are very different from 
those reported by PHOBOS (Fig.glb)). This indi- 
cates that there is apparently a significant charged 
baryon contribution in nucleus-nucleus collisions at 
forward rj region. 

Fig. shows the charg ed p ion rapidity density in 
Au-I-Au collisions RHIC and Pb+Pb collisions 
at the SPS and estimated tt" rapidity density 
from the present measurement (photon rapidity den- 
sity) at ^sjvAT = 62.4 GeV, all as a function of y- 
yheam- HIJING Calculations indicate that about 93- 
96% of photons arc from tt*' decays. From HIJING 
we obtained the ratio of the photon to 7r° yields. 
This ratio is used to estimate the tt*^ yield from the 
measured photon yield. The BRAHMS results at 
forward rapidities are slightly lower compared to the 
results from SPS energies. However, in general, the 
results show that pion production in heavy ion col- 
lisions in the fragmentation region agrees with the 
LF picture. Similar features have been observed in 
e+e^ collisions j^]. The centrality dependence of LF 
for inclusive charged hadrons and the centrality in- 
dependence of limiting fragmentation for identified 
mesons indicate that although the baryon stopping 
is different in different collision systems, the pions 
produced at forward rapidities are not affected by 
the baryon transport. 
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In summary, we have presented the first results of 
photon multiplicity measurements at RHIC in the 
pseudorapidity region 2.3 < rj < 3.7. The pseudora- 
pidity distributions of photons have been obtained 
for various centrality classes. Photon production per 
participant pair is found to be approximately inde- 
pendent of centrality in this pseudorapidity region. 
Comparison with photon and charged pion data at 
RHIC and SPS energies shows, for the first time 
in heavy ion collisions, that photons and pions fol- 
low an energy independent limiting fragmentation 
behavior, as previously found for inclusive charged 
particles. Furthermore, photons are observed to fol- 
low a centrality independent limiting fragmentation 
scenario. 
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